Implementation of three-qubit Toffoli gate in a single step 
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Single-step implementations of multi-qubit gates are generally believed to provide a simpler design, 
a faster operation, and a lower decoherence. For coupled three qubits interacting with a photon 
field, a realizable scheme for a single-step Toffoli gate is investigated. We find that the three qubit 
system can be described by four effective modified Jaynes-Cummings models in the states of two 
control qubits. Within the rotating wave approximation, the modified Jaynes-Cummings models 
are shown to be reduced to the conventional Jaynes-Cummings models with renormalized couplings 
between qubits and photon fields. A single-step Toffoli gate is shown to be realizable with tuning 
the four characteristic oscillation periods that satisfy a commensurate condition. Possible values of 
system parameters are estimated for single-step Toffli gate. From numerical calculation, further, our 
single-step Toffoli gate operation errors are discussed due to imperfections in system parameters, 
which shows that a Toffoli gate with high fidelity can be obtained by adjusting pairs of the photon- 
qubit and the qubit-qubit coupling strengthes. In addition, a decoherence effect on the Toffoli gate 
operation is discussed due to a thermal reservoir. 

PACS numbers: 03.67.Pp, 85.25.Am, 85.25.Cp 



I. INTRODUCTION 

Any quantum computing algorithm can be decom- 
posed into a sequence of single qubit rotations and entan- 
gling two-qubit gates. 1 Although this universal gate pro- 
vides a convenient and intuitive design for quantum algo- 
rithm, the design becomes so complicated as the qubits 
scale up that it may be hard to implement the quan- 
tum algorithm. The quantum gates with more than two 
qubits, thus, will be very useful for constructing quantum 
computing algorithms and quantum error correction pro- 
tocols. Physical realization of the multi-qubit gate has 
been intensively studied with various qubit systems. 

Among the multi-qubit gate, the Toffoli gate 2 
(Controlled-Controlled-NOT or C 2 -NOT gate) provides 
a universal gate set for classical computing. Recently, the 
three-qubit quantum Toffoli gate has been achieved, for 
example, with NMR qubits 3 ' 4 , ion-trap qubits, 5 and su- 
perconducting circuit £ The Toffoli gate has usually been 
implemented by decomposing it into a sequence of single 
qubit operations and CNOT gates. However, the single- 
step implementation^ of multi-qubit gate will make the 
circuit design simpler and provide shorter gate operation 
time. Further, it may achieve the lowest possible qubit 
decoherence, providing high fidelity. For the atoms in 
cavity the inter-atom interaction is mediated by the pho- 
ton mode. In this case the schemes for the single-step 
multi-qubit controlled-phase-flip gate^ and Toffoli gate^ 
have been proposed. In this study, we analyze a scheme 
for realizing the single-step implementation of Toffoli gate 
for solid-state qubits with Ising-type interaction. 

We consider a system of three qubits coupled by Ising 
interaction, where the qubit state evolution is driven by 
an external oscillating field. Two of the qubits behave as 
the control qubit, while the third one as the target. The 



external photon field is resonant with the target qubit 
for a specific control qubit state, while it is off-resonant 
for the other control qubit states. For the resonant case, 
the Hamiltonian is described by the Jaynes-Cummings 
model with a photon field, coupled vertically with the 
qubit state in a rotated coordinate. However, for the 
off-resonant case, the photon field has components par- 
allel as well as vertical to the qubit basis. In this paper, 
we analyze this modified Jaynes-Cummings model to find 
the conditions for the single-step Toffoli gate. These con- 
ditions can be obtained by investigating the commensu- 
rate oscillation periods of the target qubit states. The 
gate operation error is estimated to be so small that our 
scheme has advantages over the decomposition scheme of 
Toffoli gate. The decoherence effect of the environment 
is analyzed by introducing the interaction between the 
qubit and the thermal reservoir. 



II. MODIFIED JAYNES-CUMMINGS MODEL 

A. Effective Hamiltonian 

A qubit interacting with a microwave field, g cos tot, 
with coupling constant g can be described by a semiclas- 
sical Hamiltonian 



H, 



qubit = [-+gcosujt)a z - t q a x , 



(1) 



where e is the qubit energy splitting, t q is the tunnelling 
amplitude between different (pseudo-) spin states, and 
a z%x are the Pauli matrices. The microwave field can be 
quantized into photon field and, if we consider only the 
single-photon process, the system can be described by 



the Hamiltonian, 



H, 



qubit 



hcua/c 



e + g(a + a" 1 ") 



Tq(J x 



(2) 



where a and a^ are the photon annihilation and creation 
operators. This type of coupling can be seen, for ex- 
ample, in the circuit-QED architecture^— and the su- 
perconducting flux qubit»i2r— For three coupled qubits 
(qubits A,B,C) the Hamiltonian in the basis of {| l), | f)} 
is given by 
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Here Ei — ei + g(a + at) with i = A, B,C, and the 
Ising-type couplings are set to be equal to each other, 

jAB = jAC = jBC = J 

In Fig. [T]thc diagonal energies E ss / S » (s e {t,t}) 0I 
this Hamiltonian, excluding the qubit-photon interaction 
term g(a + at), are plotted as dotted lines. We here set 
the first two qubits as the control qubits while the third 
qubit as the target qubit. We consider the parameter 
regime, 6a,£b 3> J, so that the energy levels of coupled 
qubits | H s") are far lower while those of | tt s ") much 
higher than the energy levels of |s, —ss") as shown in Fig. 
[TJ Here, —s is opposite spin of s. Then, the tunnellings 
between the states |sss") and |s, —ss") are suppressed in 
the Hamiltonian; thus, we can set t A = and tg = 0. 
Hence, the three qubit Hamiltonian of Eq. © can be 
represented as a block-diagonal form, 



H = H\ 



tt} 



H, 



\U) 



m 



14-t) 



H, 



W>> 



(4) 



where H\, 



is explicitly given as 



haja^c 



E ss 't + E ss 'i 



Co 



E ss >f-E ss n+g(a + a^) , c 
- a z - t q a x 



(5) 



and \ss') indicates the control qubit states. 

In order to analyze the Rabi oscillation of this system 
we diagonalize the qubit part of the Hamiltonian by in- 
troducing a coordinate transformation as follows, 



V = exp [ - -a y n 



exp 



2 a y°n 



iexp 



2 a v u it 



6XP I ~~ 2° v U 



with 



tan 0j 



21° 



Ess't — E, 



^ss't 



ss' -I I 



(0) 



(7) 



Here, note that we assign j — 0,1,2,3 for \ss') = \ H 
)A It) 1 1 tt) , I tt) , respectively Then, the transformed 
Hamiltonian W — V~ 1 HV is given by 



H = U m) © U ln) © n M} © H\ U) . 



(8) 



Among the energy levels in Fig. [1] the lowest two energy 
levels are given as E^^ = —{ca + £b + ec)/2 + 3 J and 
E_n^ = — {ea + es — ec)/2 — J from the Hamiltonian in 
Eq. ©. These two energy levels are degenerate when 
ec = 4:J, and we will use this degeneracy point as the 
operating point. 

In general, the Hamiltonian is represented as a modi- 
fied Jaynes-Cummings model as follows, 



H\. 



= fouja) a ■ 
a i9 i 



1 



huj 



t)(a. 



P]9{a + 
- + *-), 



)]** 



where 



sin 0j , j3j 



flljJi 



AEl,+(2tcy 



(9) 

(10) 
(11) 



with AE SS > = E ss i^ — E ss i^, and we drop the irrelevant 
term g(a + ot)<To. This Hamiltonian has the photon field 
which has the components either vertical or parallel to 
the qubit with coupling constant a^g or fog. 

In the Hamiltonian of Eq. © the Pauli operators are 
represented in the basis of {\ss'g}, \ss'e)}, where 

\ss'g) = cos(^-/2) | ss' |) + mn.(9j/2)\a^ t), 

in(^/2)|ss' i) + cos(6j/2)\ss' t), (12) 



\ss e) 



sm 



whose energy levels are shown as solid lines in Fig. [T] 
Here, \ss'} denotes the control qubit state, and \g) and 
|e) denote the ground and excited states of the target 
qubit state. At the operating point ec = 4 J, 



AE SS , = 0, 4 J, 4 J, and 8 J 



(13) 




FIG. 1: (Color online) Dotted lines are the energy levels 
E S3 ' s n, and the solid lines correspond to the energy levels, 
(E sa , t + E 33 ,±)/2 ± (l/2)hwj, with j = 0, 1, 2, 3. The operat- 
ing point is ec = 4J. 



for \ss') = | H), | |t), I tl), an£ i | tt), respectively. Since 
the Hamiltonian Wim has degeneracy at this operating 
point, H\n) is reduced to the usual Jaynes-Cummings 
modelfii 

H lu) = hwa^a + -^-a z + |(a + a*)^ + <r_) (14) 



with 



Swo = 2ft 



(15) 



In order to describe the dynamics of qubit system, we 
introduce a rotating coordinate such as ip(t) — U(t)<f>(t), 
where 



U(t)=exp 



-II ( -Hu! (Jz 



huia* a 



(16) 



Accordingly, the Schrodinger equation Hip(i) = ih-§^ip(t) 
is written as ihM^4>{t) — Tt l (j){t) with H 1 = 
U- 1 {t)HU(t) - ihU- 1 {t){dU(t)/dt). In this interaction 
picture, the Hamiltonian T-l} „ is represented as 



ss ') 



Hi 



1 



[h(uj - uj ) + (ijg{a(t) + al(t))]a z 



+-a j9 (a(t) + a\t))(a+(t) + cr_(t)), (17) 



where the transformed operators 0(t) = U 1 (t)OU(t) 
are given by 



a±(t) =a ± e 



±iu t n f/ + \ = „Uiwt 



a 



(*) 



a' e 



f , o(t) 



(18) 



The parallel coupling term between the qubit and the 
photon field, /3jg(a(i) + cv(t))cr z , in the Hamiltonian 
%} „ can be eliminated by introducing another coor- 
dinate transformation, 



W(t) = cxp 



.J^-( a t(£)- a (*)H 

ZnujQ 



(19) 



The transformed Hamiltonian % 



\S8>) 



W- l (t)U\ ssl) W{t) - iKW- l {t){dW{t)/dt) is then 
written as 

^[ ss <> - \[K»3-^)+Pi9W)+a\t))]a* 



+ -o i s(a(t) + St(t))(er+(t) + er_(t)) 
-^•0(o(<) + «*(*)>*■ (20) 

For photon field the transformed operators 0(t) — 
W~ x O{t)W{t) are represented asi£ 

a(t) = a(t) - (0jg/2hwo)a z , 
a f (t) = <J{t) - (p jg /2fkvo)(T z , (21) 

and for spin operators^ 

er±(t) = a±(*)exp[±(& 5 //iwb)(a t (<) - a(t))]. (22) 



Then, the Hamiltonian becomes 

Wf ss /) = ^H(uj - w o)cr z + -ajg (a(t) + a f (i) 
x(cr + (*)exp[^(ot(t)-o(t))] 
+a_(t)exp[-^V(t)-a(t))]V 



hujQ 



(23) 



dropping a constant term. 

The matrix element of the Hamiltonian can be eval- 
uated in the rotating wave approximation (RWA)i£ as 
follows, [see Appendix] 

(g,n + l\ (a(t) + o f (t) - l 3 a z ) e-^ (at(t) - Q(t)) a_(i)|e,n) 



= V^+l£„(7.) + 



V" + 1 



xe-^ = F„(Ti), (24) 

(e, n| (o(t) + a f (t) - 7i o«) e^ (at(t) - a(4) V+(i)| 5 ,n + 1) 

7? 



= I V^+T£° +1 ( 7 f) + -^=[11^) + i r \( 7j 2 )] 1 

xe^EG,,^). (25) 

We can check that those two matrix elements are equiv- 
alent with each other, 



F„( 7i ) = G«(7 J ), 
by using the recurrent relations j2£ 

^( 7j 2 )-^(7l) = ^-i(7 J 2 ), 

2rl /„,2\ 



(26) 

(27) 



(n + 1)^ +1 (7|) + I'Ki^) = (n + 1)^(7/)- 

As a result, in the RWA the modified Jaynes-Cummings 
model is reduced to the conventional Jaynes-Cummings 
model with a renormalized coupling between the qubit 
and photon field. Hence, the effective Hamiltonian in 
the RWA reads 



^TT = 2 H ^ j ~ Uo)az + 2 {aa+ + a ^- ] (28) 
in the basis of {|e, n)j, \g, n + l)j}, where 

9jn = otjgFninj) (29) 

is the renormalized coupling constant. 

B. Commensurate Condition 

In this paper, we are concentrated on the ground state 

_ 1 2 

with n — 0. Then, the coupling constant cjjo = ot^ge 2 7i 
is Gaussian. Though the coupling constant obtained by 
semiclassical analysis in Ref. [2l] is oscillatory, it fits well 
with the present result for small jj because the RWA 



(?,«) 


J/t° 


g/t° 


Ql 


«3 


Pi 


ft 


7i 


73 


9io/9 


930/9 


(1,3) 


0.340 (0.358) 


0.230 (0.254) 


0.822 


0.585 


0.569 


0.811 


0.067 


0.096 


0.820 


0.582 


(2,6) 


0.460 (0.484) 


0.182 (0.199) 


0.736 


0.478 


0.677 


0.879 


0.062 


0.080 


0.734 


0.476 


(3,9) 


0.477 (0.471) 


0.128 (0.126) 


0.724 


0.464 


0.690 


0.886 


0.044 


0.057 


0.723 


0.463 



TABLE I: The values of J and g are obtained by solving Eq. (|34[) . Here the values in parenthesis are obtained from numerical 
calculation, which fit well with the analytic results for small 7 where the RWA works well, ctj(fij) corresponds to the vertical 
(parallel) coupling strength between the qubit and photon field, gjo is the renormalized coupling. 



works well for a weak coupling. For n = 0, the Hamil- 
tonian in Eq. (j2"5)l describes the quantum oscillation be- 
tween the states \g, l)j and |e, 0)^ with the frequency, 



^3 = V ("j - w o) 2 + (§30 /fy' 



(30) 



Clj is the oscillation frequency between the two states, 
1.9, l)j = \ss'g) and \e,0)j = \ss'e). 
In this study, our Hamiltonian in the RWA is written 



as 



O/ 1. RWA _ o/I.RWA 



V}I,RWA 

' H \n) 



VII, RWA 



, «;^ A ,(3D 



and we want to flip the target qubit state when the con- 
trol qubit state is | H). For j = 0, the Hamiltonian of 
Eq. (|28) is reduced to #^ A = §(acr+ + aV_), which 
describes the Rabi oscillation between the states, | 44 g) 
and | 4-4- e ) with the Rabi frequency 



fi 



R 



O = .g/fi- 



(32) 



On the other hand, for j ^ the Hamiltonian demon- 
strates a non-Rabi oscillation between the states, \ss'g) 



and 



with the oscillating frequency £L. The Toffoli 



M- 



Toffoli 



(33) 



gate requires that for a specific control qubit state the 
target qubit state flips, while for the other control qubit 
states the target qubit state remains at the original state. 
The Toffoli gate is represented as following matrix, 

10 

10 

10 

10 

10 

10 

1 

10 

in the basis of {| ft e), | ft .9>, I U e), | U 9), I it e>, \ it 
.9), I lie), 44 .9>}- 

Let us consider the initial state, |\&(0)) = (1/2) ( | 4-4- 
9) + Ut 9) + I U .9) + Itt 5))- When t = T /2 with 
T = 27r/fio, the state | 4-4 g) flips to the state | 4-4- e ), 
while the other stats \ss'g) also evolve during the time. 
Hence, in general, the Toffoli cannot be achieved at time 
t = To/2. However, if half the oscillation period, 7r/r2o, is 
an integer multiple of the other oscillation periods, 2tt/Qj 
with j y^ 0, as follows, 



H 3 



fi 2 = 2pfl , 
2qQ 



with integers p, q, the Toffoli gate can be achieved. Here, 
since A-Efj. = E^ = 4 J, ui = CJ2, cci = 0:2, and /3i =02, 
we have fii = il2- At time tk — (2k + l)7r/f2o with inte- 
ger fc, the qubit states with j = flips to the other state 
whereas the states with j — 1,2,3 come back to the orig- 
inal states. Consequently, if ilj's satisfy the conditions 
in Eq. (|3~4"|) . the Toffoli gate can be achieved at time tk- 
The coupled equations of Eq. (f3"4"|) for commensurate 
condition can be solved numerically, and partial results 
are summarized in Table [U The coupled equations de- 
termine the values of J/K and <?/ij? for a given pair of 
(p,q). Considering that the typical tunnelling frequency 
of solid-state qubits is 1-2 GHz, the interqubit coupling 
strength J/h and the Rabi frequency g/h should be of 
the order of 100MHz. The values of ctj's and /3j's are 
the prefactor of the coupling vertical and parallel to the 
qubit states, respectively, in the Hamiltonian of Eq. ©. 
<jjo is the renormalized coupling constant. We found that 
if q = 3p there are solutions for the coupled equation of 
Eq. (|3"4"]) , but we also found that there exist solutions for 
some pairs of (p, q) around the line, q = 3p. 



III. NUMERICAL ANALYSIS OF TOFFOLI 
GATE 

If the conditions for commensurate oscillation in Eq. 
(|3"4")) are satisfied, the analysis of previous section pro- 
vides a complete (maximum fidelity) Toffoli gate opera- 
tion. However, the analysis neglected the multi-photon 
processes and was performed under the RWA. In reality, 
thus, the gate operation cannot be complete. Through a 
numerical analysis with the original semi-classical Hamil- 
tonian, we can estimate the error of our Toffoli gate op- 
eration scheme. 

Similarly to the argument of section II, the three qubit 
Hamiltonian can be written as following block-diagonal 
form in a semiclassical description 



H s 



H 



ft) 



Hi 



in) 



H] 



It) 



H 



U)' 



(35) 



where H? 



F qubit 

\ss') 

rjmw 
H \ss>) 



rrqubit 



2 

g cos iota z . 



flpfj with 



t„ Or, 



E ss ,j f + E ss n 



-o- , 
(36) 



The eigenstates of H^n are given by the states \ss'g) 
(34) and \ss'e) in Eq. (1121) whose evolutions are driven by 



H! 1 ™,. Here we adjust the operating point as ec = iJ 
and set the initial states as |*(0)> = (1/2) (| LJ. g) + | J/f 
.9> + ltl.9> + ltt.9)) as before. 

The TofToli gate can be demonstrated by using the val- 
ues of J and g obtained analytically in Table fl] Here, we 
use the values for (p,q) = (1,3). In Fig. [2] (a) we can 
observe that at f2 £ = 7r the initial state | 4-4- g) evolves 
to | He), while the other states recover their original 
states as follows, 



llg) 
tig) 



lie) 

Ug) 



Its) 
tts) 



lt.9>, 

tt<?>- 



(37) 



Consequently, if the control qubit state is | 4-4-) > the target 
qubit flips between \g) and |e), while for the other control 
qubit states the target qubit remains in the original state. 
We can observe that the parameter values obtained 
by the RWA well describe the unitary evolutions for the 
Toffoli gate. The fidelity of Toffoli gate can be calculated 
by F = Tr(Af (i)M To ffoii)/8, where M{t) is the truth table 
of Toffoli gate calculated numerically. Figure[3](b) shows 
the fidelity of the Toffoli gate for (p, q) = (1, 3) and (2,6). 
We can observe that at Clot = (2fc + 1)7T the single-step 
Toffoi gate is achieved. The peak width of fidelity in Fig. 
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FIG. 2: (Color online) (a) Time evolution of the initial state 
j^(0)) by using the values of J and g obtained by the RWA 
with (p,g)=(l,3) in Table H At fit = it, the state | \\ g) 
flips to the state | 44- e )> while the other states come back to 
their initial state. As a result, the Toffoli gate is achieved in 
a single step, (b) Fidelity of the Toffoli gate for (p, q) = (l,3) 
(solid line). The fidelity approaches local maxima at Slot = 
(2k + l)ir. Dotted line shows the fidelity for (p, q) = (2,6). 




0.998 



FIG. 3: (Color online) (a) Fidelity F in the plane of (J,g) 
calculated numerically with q — 3p. The peak points are 
shown in Table [T] which corresponds to the values obtained 
analytically from Eqs. (|34[) . (b) Cut view of (a) along the 
p-axis and (c) J-axis for (p,q) = (l,3)(solid lines). Dotted 
lines show the fidelity error SF — 1 — F in log scale. 



[2] (b) is calculated as St ~ 0.5 ns at F = 0.999 with 
iJ?=lGHz, which is sufficiently wide for current state-of- 
the-art technology. 

Figure[3] (a) shows the fidelity F in the plane of (J, g). 
The peak points denote the local maxima of fidelity for 
q = 3p, which are listed in Table fl] Since the RWA works 
well for weak coupling, the RWA results fit well with 
the numerical calculation for small g in Table [I] Figure 
13] (b) and (c) show the cut views of fidelity in (a) with 
(p,q) — (1,3) along the g-axis and J-axis, respectively, 
which pass the local maximum. The peak width with 
F = 0.999 is Sg/tf « 0.008 and 6J/t° ~ 0.01. For 



c _ 



1 - 2GHz, Sg and 8 J are of the order of 10MHz. 



t 

In Fig. [3] (b) and (c) the fidelity errors 5F = 1 - F are 
shown. The maximum fidelity (F=l) can be achieved by 
using the Jaynes-Cummings Hamiltonian of Eq. (|28|l in 
the RWA. In fact, however, our numerical calculation is 
based on the semiclassical Hamiltonian of Eq. (|T|), which 
includes multi-photon processes and fast mode oscilla- 
tions. As a result, due to these imperfections we have 
the fidelity error; SF ~3x 10 -4 which is 3 times larger 
than SF ~ 10~ 4 for the two-qubit oscillations^! The 3- 
qubit Toffoli gate can be implemented by decomposing 
it into a sequence of 10 singlc-qubit operations and 6 
controlled-NOT operations.— Although the fidelity error 
for the present three-qubit oscillation is larger compared 
to two-qubit gate operation, our single-step implementa- 
tion scheme has an advantage in this respect. 



IV. DECOHERENCE ANALYSIS 



The environment surrounding the qubit system is 
known to invoke a decoherence in the qubit state. The 
environment can be described as a reservoir of thermal 
bath. We consider that our qubit of Eq. (|28)l is coupled 
with the reservoir. If we trace out the reservoir degree of 
freedom from the total density matrix, the time evolution 
of the reduced density matrix p is given by the Master 
equation in the Born-Markov approximation^ - — 



dp i r^LRWA 1 



7 r^I.RWA r^I.RWA 11 / 00 \ 

^ U \ss-) >i n \ss>) .P]]. ( 38 ) 



where 7 depends on the temperature and the spectral 
density of the reservoir. 

Here, if we introduce a dressed state basis, 



COS( 



J" 



\9,n- 



1) 



sin< 



]n\ 



\-n) = -sm<f> jn \g,n+ 1} + 



cost 



jn\ 



with 



tanc 



9jnVn + 1 



(39) 



=,(40) 



H(u>j - cj )+^gj n (n + l) + h 2 (uj j - w ) 2 

can be diagonalized 



the qubit Hamiltonian of Eq 

as 



U]?^\±n)=Ef n \±) 



(41) 



Ef n = fiw (n + ± ly/gKn + V + Vfa-uo) 2 . 

In this dressed state basis the Master equation of Eq. 
jSJ) can be exactly solvable, providing the time evolution 
of the density matrix — 



rPP 



dV 



(*) = 



^(^-<)+^(^-<) : 



pZ.it), 



(42) 



where p,j/ G { 
P 



-} and 



Pnm I + n) (+m\ + p£ m I + n) (-m\ 

+Pn™\ -n)(+m\ +p~~\ -n)(-m\. (43) 



Let us consider that initially the target qubit is at the 
ground state, \ip(0)) — \g,n+l). The density matrix 
p(0) = \g, n+l){g, n+1 | can be represented in the dressed 
state basis by using Eq. (I39|) . Then, the density matrix 
of initial state is given by 



Pnn(0)=COS (f> jn , p+„(0) 
PnnW = -Sin^„ 



,J jn, Hnn 



(0) 



(44) 



By tracing out the density matrix over the photon field 
we obtain the density matrix of the qubit system T = 
^2i(l\p\l), where T is a density matrix in the basis of 




FIG. 4: (Color online) Fidelity of Toffoli gate in the presence 
of the reservoir for several qubit-reservoir coupling 7's. 



{|e), \g)}. The expectation value of ground state popula- 
tion Qss'g(t) = — Tr[Tc z ] is calculated as 

Qss' S (t) = cos 2 2<j>j n + sin 2 2<f>j n cos(fl,„£)e~ 7f ^™*, (45) 

where H jn = (E+ n - Ej n )/h. 

In our scheme, for the control qubit state | 44) (j = 
0), the target qubit states are resonant with the photon 
field. For the resonant case the population expectation 
has been derived^ In our case Qn g (t) is reduced to 

Q Ug (t) = cos(g „V^TTt/^)e-^- ( " +1)t/?l2 , (46) 

which is consistent with the previous result. 

Since we are considering the case n = in this study, 
we have 

Q ss , s (t) = cos 2 2<f>j + sin 2 2<j>j cos(rj J i)e _7 °? t (47) 

with ilj = Qjo and <pj = 4>jq. Here, we introduce the 
occupation probabilities 

Pss> g (t)= 1+Qs 2 S ' S ^ and P ss , e (t) = l-P ss , g (t). (48) 

For 7 = these probabilities correspond to the numerical 
result of Fig. [5] (a), where the initial target qubit state 
is the ground state. 

The fidelity for the Toffoli gate is given by F = 
Tr(M(i)M To ffoii)/8, where the truth table matrix M(t) 
has a block diagonal form 

M{t) = M(% t) 8 M{t)\ U) 9 M(t) m) © M(t) lu) (49) 

with 

M(t)\ A = ' *\ss'e)^\ss'e){t) " \bs' e)-¥\sa' g) \f) 



P, 



ss' g)— >\ss'e) 



(*) P\ 



\as' g)-*\es' g) 



(*) 



(50) 



Here, if initially the qubit evolves from the ground 
(excited) state, P\ SS 'g(e))->\sa'g(e))(t) i s the probabil- 
ity that the qubit will occupy the ground (ex- 
cited) state at time t. Further, we have the rela- 
tions that P\ ss 'g)^\ss'g)(t) = P\ s8 'e)^\ss'e)(t) = 1 ~ 

P\ss' g )^\ss'e)(t) = 1 - P] ss >e)->\ S 8>g) (*)• The probability 
P\ss'g(e))->-\ S s'g(e))(t) is given by the occupation probabil- 
ity of Eq. fltSJ} as 



P\ss'g)->\ss'g)(t) — P\ss'e)^\sa'e) (*) — 5 (51) 

and, thus, the fidelity is written by 

^ = ~[4 + QnsW + <W*) + <9;tg(*) - Qus(*)l- (52) 

We show the damping of fidelity in Fig. @]for several 7's. 
As shown in Fig. |4l the maximum fidelity decays with 
the operation time. 



V. SUMMARY 

We study a scheme for the single-step Toffoli gate for 
three qubits coupled by the Ising interaction. The photon 
field is resonant with target qubit for a specific control 
qubit state whereas it is off-resonant for the other con- 
trol qubit states. For the target qubit state with a reso- 
nant photon field, the qubit-photon system is described 
by the Jaynes-Cummings model demonstrating a Rabi 
oscillation, while for the other control qubit states it is 
described by a modified Jaynes-Cummings model. We 
found that in the RWA the modified Jaynes-Cummings 
model is reduced to the usual Jaynes-Cummings model 
with a renormalized coupling between qubit and photon. 
The single-step Toffoli gate can be achieved for a 7r/2 
rotation, if the oscillation periods of target qubit satisfy 
commensurate conditions. The commensurate condition 
determines the values of the qubit-photon coupling g and 
the qubit-qubit coupling J for achieving the Toffoli gate. 
These values fit well with those obtained by numerical 
calculation for weak coupling g . The fidelity is shown 
to be high and the peak width of the fidelity to be wide 
enough for implementing the Toffoli gate. The fidelity er- 
ror is shown to be small, and thus our scheme has advan- 
tages over the decomposition scheme. The decoherence 
effect from the environment is discussed. 
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Appendix A 

According to the Baker-Campbell-Hausdorff formula 
one of the terms in Hamiltonian of Eq. (|23[) . for example, 



a(t)e 



- 7j -(a+(i)-a(t)) 



a_(t) 



(Al) 



E (- 7j aW( 73 . a p e _ 

' J r\ el 



where 



7j 



_ yj 



Pj9 



hbJ 



(A2) 



In this case, by using Eq. (|18[) the time dependent part 
of the term is given by e Hr-s-i)wt-iu> t ^ In the RWA fast 

oscillating terms are neglected. Hence, in order for the 
term to survive in the RWA when uj = ojq, it should be 
that r — s = 2. This means that for this term we need 
to evaluate (n + 2\e~ Jj ( a (*)~ a (*))|n) and then the only 
non- vanishing matrix element is 

( ff ,n + l|a(i)e^ (at(t) - a(4) V-_(t)|e,n). (A3) 

Here, \g) and |e) denote the qubit states, and \n) and 
I n + 1) the photon number states. For another term 

a t (i)e~T^ at (*)~ Q ( t )) C r_(i), it should be that r = s, and 

we need to evaluate (n\e~ lj ( a w -a W)|n). Then, the only 
non- vanishing element is 

{ff,n+l|o t (t)c^ (otw - a(t) V_(t)|e,n). (A4) 

The terms with <r+ (t) and a z in the Hamiltonian of Eq. 
(|23[) can also be evaluated in a similar manner. 

These matrix elements can be evaluated by the formula 
for the displaced number stated for n < m, 



{n\e iai -^ a \m) = (-7*)™""-! —^L^-^d^e^ 1 ^ 2 , 



(A5) 



m |e-T ot +T*°|n) 



and for n > m, (n\e la 7 a \m) 
where L^(x) is the associated Laguerre polynomial. 
By summing up all the contributions we have 

(g,n+l\ (o(t) + a+(t) - 7j a z ) e"'W< at W- < t »CT_(i)|e,n} 
= L^TlL^ 2 ) + -0—[LWI) - Li(tf)]\ 

x e~^' 2 = F„( 7i ), (A6) 

(e, n\ (a(t) + a\t) - l 3 a z ) e^^^' a ^a + (t)\g,n + 1) 



I V^TTi° +1 ( 7j 2 ) + -/=4=[l;U(7|) + Ll(i 2 )] I 

xe^ 2 / 2 ^G„( 7j ). (A7) 
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